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O ■ Abstract 



We study the solar neutrino problem within the framework of a parametrized 
post-Newtonian formulation for the gravitational interaction of the neutrinos, 
which incorporates a violation to the equivalence principle (VEP). Using the 
current data on the rates and the energy spectrum we find two possible os- 
cillation solutions, both for a large mixing angle. One of them involves the 

X' 

$—1 ' MSW effect in matter and the other corresponds to vacuum oscillations. An 

interesting characteristic of this mechanism is that it predicts a semi-annual 

variation of the neutrino flux. Our analysis provides new constraints for some 

VEP parameters. 
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I. INTRODUCTION 

Several experiments sensitive to solar neutrinos have measured a v e flux, with results 
lower than the values predicted by standard solar models (SSM) for different neutrino ener- 
gies: the Homestake CI radiochemical experiment , with sensitivity down to the lower part 
of the 8 B spectrum and the 7 Be line, the two radiochemical 71 Ga experiments, GALLEX |2j 
and SAGE [^], which are sensitive to the low energy pp neutrinos and above, and the water 
Cerenkov experiments, Kamiokande Q and SuperKamiokande (SK) [EJ, which can observe 
only the highest energy 8 B neutrinos. A combination of any two of the experiments disfavors 
an astrophysical solution to the problem, and seems to indicate that a non-standard physical 
process is modifying the energy spectrum of solar neutrinos. 

A widely accepted explanation of the discrepancy is based on the assumption that non- 
degenerate massive neutrinos do undergo flavor oscillations, either in vacuum or within the 
Sun (MSW effect) ||. Another less orthodox mechanism for neutrino oscillations, which 
does not need neutrinos to have a mass, was proposed several years ago |7j and requires 
the coupling of neutrinos to gravity to be flavor dependent, i.e., a violation of the equiva- 
lence principle (VEP) in the neutrino sector. Some phenomenological consequences of this 
mechanism have been examined in a number of papers p|-p^| . 

In a recent work |D|] we developed a generalized VEP mechanism for neutrino oscillations, 
which is based on an extended parametrized post-Newtonian formalism (PPN). Here we 
apply this approach to the concrete situation of solar neutrinos, and in particular to the 
analysis of the seasonal variation of the signal. Using the latest data on total rates from the 
five experiments, and those on the energy spectrum and the seasonal variations from SK, 
we determine the allowed regions and the best-fits values for the oscillation parameters. We 
show that a solution to the solar neutrino problem is possible within the VEP scheme, not 
only for MSW matter-enhanced transformations but also for vacuum oscillations. 

In the solar system the gravitational field receives contributions from several sources. 
Assuming, as is commonly done, that the potential vanishes at an infinite distance from the 



source, the dominant contribution is given by the Great Attractor, with small perturbations 
produced by galactic clusters, our galaxy, and the Sun. Consequently, it is reasonable to 
approximate the potential by a constant of the order of 10~ 5 |l6|j . The effect of this potential 
regarding a possible VEP mechanism has already been analyzed by Halprin, Leung, and 
Pantaleone |Tj|. Here we follow the more general approach of Ref. jT5|], which incorporates 



not only a possible flavor dependence of the gravitational couplings, but also the most general 
violations to Einstein gravity in the context of metric theories. In this approach, the metric 
is given by the Minkowskian one plus source dependent perturbations: g^ u = r/^ + h^ v . 

The assumptions for constructing the metric in the PPN formalism involve virialized 
sources such that ^ ~ w 2 , where the quantities M, R, and w represent estimations of the 
order of magnitude of the mass, distance and characteristic velocity of the source. In what 
follows we keep only first order corrections to the flat space-time metric r]^, and we neglect 
a possible angular momentum of the Great Attractor, which in any case would lead to very 
small corrections. Thus we have h oi = 0, while the non-null corrections are given by 

h 00 = 2iU + 0(w A ), (1) 

h i:j = 2-fUS i:j + TUij + Oiw 4 ) , (2) 

where 7, 7', and V are adimensional parameters of the PPN expansion (up to order w 3 ). 

In the particular case of Einstein gravity we have T = 0, and 7 = 7' = 1 . The potentials 

responsible for the metric perturbations are 

u= f p(r') d 3 r' JJ = r p(r')(r^ - rQfc - r\) d 3 r' 

J I r — r' I lJ J I r — r' | 3 

with p(r) being the mass density of the source of the gravitational field. We are using a 

system of unities with G — h — c — 1. 

For a confined and distant source, U can be approximated by 

If we take the z-axis along the direction determined by the solar system and the gravitational 
source, we then have U zz ~ U. However, the components U xz and U yz are proportional to 



A9U, where A9 is the angular size of the source, while U xx , U yy , and U xy are of the order 
of (A0) U. Since the Great Attractor is a rather extended object with an angular size 
of the order of 10 _1 |[L7|, in the case of the solar system there are only three relevant 



types of Uij contributions: (i) those coming from our galaxy, which are of order 10~ 6 , 
(ii) a longitudinal component from the Great Attractor, of order U zz ~ U ~ 10 -5 , and 
(hi) transverse-longitudinal components also produced by the Great Attractor, of the same 
order as the galactic contributions, U xz ~ U yz ~ 1CT 6 . Therefore, possible VEP flavor 
oscillations of solar neutrinos would be characterized by three main effects: an isotropic 
effect (U ~ 10~ 5 ), and two anisotropic effects (U zz ~ 10~ 5 , Uij ~ 10~ 6 ). In the next section 
we review the essential ingredients of the VEP mechanism for neutrino oscillations within 
the context of the PPN formalism, and in Section III we apply it to the study of the solar 
neutrino problem. 

II. VEP INDUCED OSCILLATIONS 

For simplicity, in what follows we consider that there are only two (massless) neutrino 
flavors, v e and v^. In our VEP scenario they are assumed to be linear superpositions of the 
gravitational eigenstates z/f and v^ , with a mixing angle 9 g . Each gravitational eigenstate 
is characterized by a different set of PPN parameters, {7 a ,7 /a ,r a } (a = 1,2). This leads to 
different dispersion relations for the i/f , which can be approximated by [|l5j 



E a = p 



1 _^ + 1 ap_ r a u .m 



p 2 



(5) 



Suppose that the initial state produced at time to corresponds to a pure electron neutrino. 
Then, for a constant gravitational field the survival probability after traveling a distance 
L = t — t is 

P{v e -> u e ) = 1 - sin 2 26 g sin 2 — . (6) 

According to this, neutrino oscillations will appear whenever a non null mixing is generated 



because of flavor dependent gravitational interactions. These oscillations have a character- 
istic length given by 

A, = £7, (7) 

l A o| 

with 

PiPj 



Aq = -C/2 — E\ = — E 



(Sj' + 8^)U + STU k 



(8) 
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where E ~ p is the neutrino beam energy, and 

5 7 = 7 2 - 7 1 , Sj'^j^-j' 1 , 5T = T 2 -T 1 . (9) 

In contrast to the ordinary vacuum oscillations induced by a mass difference, where X m = 
AirE/Sm 2 is proportional to the energy, the effect we are considering here has an oscillation 
length that goes with E~ l . This leads to observable distinctions between both mechanisms 
and makes the gravitational induced oscillations suitable to be observed with higher energy 
neutrinos P,|T2||. Note that even though the overall sign of the gravitational potential is 



irrelevant for oscillations, the relative signs among differences of the PPN parameters are 
very significant. If we assume that these differences are all of the same order, then the most 
important directional effect would be given by the quadrupolar contribution corresponding 
to U zz . 

As is well known, flavor transformations of massive neutrinos are affected by their in- 



teractions with matter | lq| . Neutral current interactions are flavor diagonal and can be 



ignored, as long as we do not consider sterile neutrinos and neutrinos are not part of the 



medium |JL9|| , but this is not true for the charged current interactions. As a consequence, the 
forward scattering amplitude is not flavor diagonal and depends on the leptonic content of 
the matter, which gives place to important consequences such as the MSW effect. A similar 
phenomenon happens for the VEP mechanism in the presence of matter. In this case, the 
flavor evolution for relativistic neutrinos propagating through a constant gravitational field 
is governed by the equation 



a v > 



i— 
dt 



/„ \ 



H(t) 



VV 



(10) 
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where the Hamiltonian 7i(t), after discarding an irrelevant overall phase, can be written as 



Hit) 



Ao 
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-cos29 g sin20 g ' 



V 



+ 



b(t) 



sm28g cos29 g i 
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(id 



The first term arises from VEP and the second term accounts for the matter effects on 
the neutrino propagation. For a normal matter background, as in the case of the Sun, we 
have b(t) = v2GViV e (£), where Gf is the Fermi constant and N e (t) denotes the electron 
number density. An extended gravity like the one here considered could also affect the 
electroweak Lagrangian, but the combined effect should be of the order UGf- Therefore, 
they are strongly suppressed and we do not include them in our discussion. 

Taking into account the dominant contributions due to U and U zz ~ U, the coefficient 
A reduces to 



A r; 



-EU 



O7 + <V) + <5r cos 2 D cos 2 (a - A) 



(12) 



Here, a is the right ascension of the Sun, and A and D are the right ascension and decli- 
nation of the Great Attractor in ecliptic coordinates. The second term in A arises from 
the quadrupolar potential of the gravitational source and generates a seasonal dependence 
in the oscillation wavelength, as first discussed in Ref. |l]|. This effect went unnoticed in 
previous work on the subject [|9|-|T3ll, where only the contribution coming from the Newto- 
nian gravitational potential was considered. To isolate the anisotropic contribution, it is 
convenient to reparametrize A as follows 



A = -EU5^i 



1 + cos 2 {a- A)--) 8 



(13) 



where £7 = (£7 + Sj')/(1 — 5/2) and S = 5T cos 2 D/S'y, so that the annual average of Ao is 
independent of 8. We will define 5 positive, because (8, A) is equivalent to (—5, A + 7r/2). 

At any time, Ti(t) can be diagonalized by a unitary transformation characterized by an 
angle 9 m (t) 
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sin20 m (*) = A ° Sin2gg • (14) 

^(A cosW g - 6 e (t)) 2 + (A sin2^) 2 

There exists a resonant flavor conversion when the diagonal elements of the Hamiltonian 
vanish, i.e., when 

V2G F N e (t R ) = A cos 26 g , (15) 

and in this case the mixing in matter is maximal (sin2# m = 1). 

The efficiency of the conversion mechanism depends on the adiabaticity of the process. 
For a constant gravitational field, the average probability for a u e produced in the Sun to 
reach the Earth reads 

P(y e - ^ e ) = i + 1(1 - 2P c )cos20^cos20 fl , (16) 

with (P m = 9 m (to). The function P c represents the probability of transition between the 
instantaneous eigenstates of Ti(t). It embodies the total correction to the adiabatic result 
for (P Ue ), which corresponds to P c = 0. 

Except for regions close to the center and the surface, the electron density in the Sun is 
well approximated by an exponential profile |2(| . Thus the change of the electron density 



along the path of a neutrino moving radially within the Sun can be written as 

N e (t) = N e (t )exp[-(t-t )/r }, t > t (17) 

where N e (t ) is the density at the production point and r is a parameter to be adjusted 
according to the region |21|| . In the SSM N e takes its maximal value at the center of the 
Sun, where it is approximately equal to 100 Na g/cm 3 , where Na is the Avogadro number. 
For N e (t) as given in Eq. (|17D , the following formula for P c has been derived in a given 



approximation from the exact analytical solution of the evolution equation p2 

/ cos26> g \ 
V l-cos 26» 9 J 



cxp 



TXK 



exp 
where the adiabatic parameter k is 



/ 2cos2fl„ \ 
1XK \ sin^ 29 g ) 



(18) 



|A„| Sin2 ^ tan2 5 (19) 



Ne(t R ) 



dN e (t) 
dt 



tn 



In the denominator of the last formula we have discarded any term associated with variations 
of the gravitational field with the distance. From the above expressions, we see that P{y e -^ 
u e ) depends on the electron density in the production zone and the logarithmic derivative 
of the density in the transition layer. 

For k ^> 1, P c is exponentially small. On the other hand, when k < 1 there are 
considerable corrections to the adiabatic approximation that reduces the magnitude of the 
resonant transformation. Nonadiabatic effects become important when k is of order 1, 
provided that the neutrinos go through a resonance. If b(to) < Aocos2#, level crossing cannot 
occur, P c = and neutrino propagation will be adiabatic even for k < 1. An effective way 
to account for this situation is to multiply the expression of Eq. ([T^) by the step function 
O(6(to) — Aocos2#), so that the transition probability vanishes when neutrinos are produced 
below the resonance. The MSW survival probability (P Ue ), averaged on the production 
region for the 8 B neutrinos, is plotted in Fig. 1 as a function of Ao for different mixing 
angles. These curves have been obtained using the electron density predicted by the SSM 
20| , |23[| . Notice that, in contrast with the common MSW mechanism for massive neutrinos, 



in the case of VEP the adiabatic edge is at higher energies, whereas the nonadiabatic edge 



is at lower energies [12|. The adiabatic edge shifts towards higher energies with decreasing 



mixing angle, as seen from Eq. (|T9"D and the condition k ^> 1. 

III. NEUTRINO EVENT RATE AND ENERGY SPECTRUM 

In the presence of neutrino oscillations, the capture rate for the radiochemical experi- 
ments, such as 37 C1 and 71 Ga, is given by: 

R(E) = g(t) ]T / dE v $ k (E v ) (P Ue ) o{E v )dE v , (20) 

k Jo 

where a(E u ) is the cross section for neutrino capture and <&k{E v ) is the fc-component of 
neutrino flux spectrum. Here, g(t) is a geometrical factor due to the Earth's orbit eccentricity 
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and (P Ue ) is the survival probability averaged over the production regions for the different 
neutrino sources. 

For neutrino-electron scattering experiments, such as SK, the solar neutrino induced 
event rate can be written: 



/OO ("OO 

dE e E(E e ,E) / dE v $(E, 



da e (E u ,E e ) _ da^(E u , E e ) 



(21) 



dE e dE e 

where E v is the energy of the incident neutrino, E e is the electron kinetic energy, and <&(E V ) 
gives the neutrino flux spectrum. The function S(E e , E) characterizes the Superkamiokande 
efficiency to measure the energy of the scattered electrons [||, and dae/dE e (£ = e,/z) is the 
differential cross section for the vg — e elastic scattering, where E e is the electron kinetic 
energy. This differential cross section can be calculated from the electroweak theory, and is 
given by 

dot cr 



dE P m P 



J e 



2 . 2 (-, E e\ 2 (m e E e S 



(22) 



with cr = 8.8 x lCT 45 cm 2 , gn = sin 2 9w, and gi = ±| + sin 2 6 W . The upper sign corresponds 
to v e — e and the lower sign to v^ — e scattering, respectively. For the energy interval of 
solar neutrinos da^/dEg = (0.155 — 0.166)da e /dE e . 

The VEP mechanism begins to be significant when half of an oscillation is about equal 
to the Sun-Earth distance. According to Eqs. (0) and (0), for a 10 MeV neutrino this 
corresponds to \USj\ ~ 10~ 25 , in which case we have pure vacuum oscillations. For larger 
values of \USj\ the oscillation wavelength shortens, and when it becomes smaller than the 
solar radius the effect of the background matter turns out to be relevant through the MSW 
effect, with the mixing angle 8 m given by Eq. (|I4|) . To compute the event rate we follow in 
general the scheme of Ref . [0] . The ingredients used in our computation have been developed 
in different places. The matter effects on the calculation of (P Ue ) were incorporated by 
applying the analytic formula given by Eqs. ([TBI) and (|l8l), as discussed in Ref. plj . The 
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electron density is given in Refs. [EOl and p3|, while the cross sections and the neutrino 



fluxes were obtained from Refs. |20| and 

We identify three regions in the \U5 : y\-sm' 2 29 parameter space for the VEP induced 
oscillations which are compatible with the observed total rates. Two of them correspond 
to MSW-enhanced VEP oscillations, whereas the third one is associated to vacuum VEP 
oscillations. The MSW solutions and the vacuum oscillation solution are separated by three 
orders of magnitude in |£/#7|. To identify these regions we use a standard x 2 analysis f25| of 
the data from all the solar neutrino experiments, taking into account both the experimental 
and theoretical errors. 

As Fig. 2 shows, for the MSW VEP oscillations there are two 99% c.l. regions allowed 
by the measured rates in all the experiments. One of them is a small mixing angle solution, 
with 3.2 x 1(T 3 < sin 2 (2# 9 ) < 5.7 x 1(T 3 and \USj\ ~ 3.2 x KT 19 , and the other is a large 
mixing angle solution, with 0.6 < sin 2 (26» g ) < 1 and 1CT 22 < \USj\ <4x KT 21 . The best 
fit for the small mixing angle is obtained with sin 2 (26' 9 ) = 4 x 10~ 3 , whereas in the case of 
the large mixing angle it occurs at \USj\ = 1.58 x 10~ 22 and sm 2 (29 g ) = 0.87. At 94 c.l. the 
small mixing angle region disappears, and only the large mixing angle region remains. 

Our analysis reveals that there is another allowed region, which corresponds to vacuum 
VEP oscillations and is shown in Fig. 3. At 99% c.l. the main sector is bounded by 
0.75 < sin 2 (2# 9 ) < 1 and 10~ 24 < \U6j\ < 10~ 22 . The values of the parameters for the 
best-fit point are \U6j\ = 1.82 x 10~ 24 and sin 2 (2# 9 ) = 1. The MSW VEP solutions are 
consistent with those already found using the Newtonian approximation for the gravitational 
interaction fl2]]l3|| , while the new solution given by the vacuum VEP oscillations has been 
independently derived in a recent work [p6fl . In previous studies it has been argued that 
when half of an oscillation corresponds to the Sun-Earth distance for 10 MeV neutrinos 
the 8 B neutrinos are depleted but the lower-energy 7 Be neutrinos are unaffected, which is in 
contradiction with the experimental data P,|T2"|| . However, a good agreement can be obtained 



if the wavelength is tuned for an energy E t close to the energy of the Be line. In this way we 
have the required suppression of the lower-energy neutrinos, and due to the inverse energy 
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dependence of the oscillation length \ g for VEP oscillations, we can also have a reduction 
by about 50% of the neutrino flux for higher energies, E^ = nE t with n integer. This is the 
origin of our vacuum VEP solution, with \ g tuned for 1.13 MeV neutrinos. 

Besides the total rate, the SK collaboration has provided spectral information on the 8 B 
solar neutrinos [27|. The measured energy spectrum of the scattered electrons is divided into 



bins having a width of 0.5 MeV in the range from 5.5 MeV to 14 MeV. An additional bin 
comprises of events with energy from 14 MeV to 20 MeV. Figs. 4 and 5 show the x 2 analysis 
for the energy spectrum || corresponding to the MSW and vacuum VEP oscillations, and 
in Fig. 6 we display the spectrum of the best VEP solutions together with the experimental 
data. The small-angle MSW solution is excluded by the energy spectrum at 99% c.L, 
while both the vacuum solution and the large-angle MSW solution are allowed at 90% 
c.L Figs. 7 and 8 display the x 2 analysis carried out with the whole set of data, including 
simultaneously the total rate and the SK spectrum measurements, with the individual x 2 
treated as independent 0. 

The eccentricity of the Earth's orbit produces a geometrical 7% variation of the neutrino 
flux since the Earth-Sun distance changes throughout the year. Due to the dependence of 
(P Ue ) on distance, an anomalous additional effect can be caused by the presence of the usual 
vacuum oscillations between massive neutrinos. Both effects are characterized by a one 
year period. Some indications of a seasonal variation in the neutrino flux from the Sun has 



already been seen in the GALLEX and Homestake experiments |28j . In Ref . [£7| , SK has also 
presented preliminary results that slightly favor a seasonal variation of the solar neutrino 
flux for E e > 11.5 MeV in addition to the geometric variation. Within the present VEP 
oscillation scheme a non-geometrical seasonal variation of the flux is caused by the presence 
of the term proportional to 5T in A (see Eq. (|12"D), which would produce a six month period 



variation. As a consequence, in contrast with the usual mass mechanism, the effect should 



be observed even in the case of MSW transformations. The authors of Ref. |26|] conclude 
that no strong seasonal variation in the solar neutrino signal is expected for vacuum VEP 
oscillations. The difference with our result is due to the fact that in their analysis they follow 
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the common prescription to incorporate gravitational effects only through the Newtonian 
potential. 

The seasonal variations of the flux above 11.5 MeV predicted by the best-fit VEP solu- 
tions including the anisotropic term are shown in Fig. 9, together with the SK data. All 
the solutions have a similar behavior and are compatible with the data within the present 
statistical accuracy. This is a consequence of the time resolution of the present experimental 
results. In principle, these solutions could be easily discriminated if the time resolution was 
improved, because their actual temporal dependence is very different, as Fig. 10 shows. The 
X 2 analysis based on the data of SK for energies above 6.5 MeV and 11.5 MeV are shown 
in Fig. 11, in terms of the parameters 5 A and S, where 5 A denotes the difference between 
the perihelion right ascension of the Sun and A modulo it. The best-fit solutions are: (a) 
5 A ~ 140° for MSW oscillations with large mixing angle, (b) S A ~ 60° for MSW oscillations 
with small mixing angle, and (c) 5 A ~ 50° for vacuum oscillations. Since 30° and 150° are 
the values of 5 A that are consistent with the position of the Great Attractor, the previous 
results seem to favour the large-angle MSW solution. This is not very conclusive because of 
the poor angular resolution of the data and the uncertainty concerning the position of the 
Great Attractor. Nevertheless, the analysis does give an improved boundary for the possible 
values of the parameter 5. At a 90% c.l., we have 5 < 0.09 for the small angle solution, and 
5 < 2 for the large angle solution. 

IV. FINAL REMARKS 

We have used the present experimental result to re-examine the possibility that the 
VEP mechanism can provide a consistent solution to the solar neutrino problem. The total 
neutrino rates give three allowed regions: two of them correspond to a MSW solution, and 
the third one to a vacuum solution. The small-angle MSW solution is excluded at 94% c.l. 
and the large-angle MSW solution is discarded at 88% c.l. The most favored solution is 
given by the long- wavelength vacuum oscillations, whose best fit has a 35% c.l. The three 
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VEP solutions predict a seasonal dependence of the neutrino flux, in good agreement with 
the data. An improved time resolution is necessary in order to discriminate between the 
different periodic behaviors. The most restrictive conditions arise from the energy spectrum. 
In this case, the small-angle MSW solution is clearly ruled out by the SK data whereas the 
vacuum VEP oscillations are favored. The MSW oscillations with large mixing angle also 
remain as a possible solution. 

More accurate data are necessary to properly establish the viability of the VEP mecha- 
nism as an adequate explanation of the solar neutrino deficit. However, the present analysis 
is sufficient to set new boundaries on the PPN coefficients that parametrize the violation of 
the equivalence principle. An improved resolution in the spectral measurements is required 
to establish the existence of a six-month period variation in the 8 B neutrino flux, which is a 
signature of the VEP mechanism and makes a clear difference with other possible solutions, 



such as the standard vacuum oscillations of massive neutrinos p9 |. 

Taking into account the different energy dependence of the VEP and mass mechanisms, 
a combination of both oscillations would give an excellent agreement with the experimental 
data. This possibility will be explored in a forthcoming paper. 
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FIGURE CAPTIONS 

Fig. 1: Neutrino survival probability for MSW VEP oscillations as a function of A for 
different mixing angles. The probability has been averaged over the production region of 
the 8 B neutrino. 

Fig. 2: MSW solutions, total rates only. Allowed regions at 99% c.l. and 90% c.l. 
(darker shaded region) in the U5'y-sm 2 (29g) parameter space. The best-fit point is indicated 
by a cross. 

Fig. 3: Vacuum oscillations, total rates only. Allowed regions at 99% c.l. and 90% c.l. 
(darker shaded region) in the U5^f-sm 2 (29 g ) parameter space. The best-fit point is indicated 
by a cross. 

Fig. 4: MSW solutions, SK spectrum only. Excluded regions at 99% c.l. (darker shaded 
region) and 90% c.l. in the U5 : y-sm 2 (29 g ) parameter space. The best-fit point is indicated 
by a cross. 

Fig. 5: Vacuum oscillations, SK spectrum only. Excluded regions at 99% c.l. (darker 
shaded region) and 90% c.l. in the f/57-sin 2 (2# 9 ) parameter space. The best-fit point is 
indicated by a cross. 

Fig. 6: SK measured energy spectrum and best fits for the small mixing angle (SMA) 
MSW solution, large mixing angle (LMA) MSW solution, and vacuum oscillation (VO) 
solution. The solid line corresponds to the VO solution for the best fit of the spectrum data 
and the long-dashed line is the VO solution for the combined best fit of both the total rates 
and the spectrum data. 

Fig. 7: MSW solution, total rates and energy spectrum. Allowed regions at 99% c.l. and 
90% c.l. (darker shaded region) in the U5j-sm 2 (29 g ) parameter space. 

Fig. 8: Vacuum oscillations, total rates and energy spectrum. Allowed regions at 99% 
c.l. and 90% c.l. (darker shaded region) in the U5 : y-sm 2 (29 g ) parameter space. The best-fit 
point is indicated by a cross. 

Fig. 9: Seasonal variation of the flux above 11.5 MeV. We have plotted eight annual 
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bins for the different best-fit VEP solutions and for the SK data (black squares). The 
anisotropic effects due to ST have been considered in the calculation. The solid line shows 
the geometrical variation due to the eccentricity of the Earth's orbit. 

Fig. 10: Actual temporal variations of the flux for the different best-fit VEP solutions 
compared with the SK data. 

Fig.ll: Allowed 90% c.l and 35% c.l regions in the (5-8 A) plane for anisotropic VEP 
solutions: (a) VO solution, (b) SMA MSW solution, and (c) LMA MSW solution. The 
best-fit points are indicated by crosses. 
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